Nonlocal Excitations and 1/8 Singularity in Cuprates 
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Momentum-dependent excitation spectra of the two-dimensional Hubbard model on the square 
lattice have been investigated at zero temperature on the basis of the full self-consistent projection 
operator method in order to clarify nonlocal effects of electron correlations on the spectra. It is found 
that intersite antiferromagnetic correlations cause shadow bands and enhance the Mott-Hubbard 
splittings near the half-filling. Furthermore nonlocal excitations are shown to move the critical 
doping concentration 5^, at which the singular quasiparticle peak is located just on the Fermi level, 
from <5j* = 0.153 (the single-site value) to 5^ = 0.123. The latter suggests the occurance of an 
instability such as the stripe at 5% = 1/8. 
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I. INTRODUCTION 

Low energy excitations and associated unusual behav- 
iors of electrons in two-dimensional copper-oxide super- 
conductors have been a question under debates in the 
past quarter century. Although theoretical studies based 
on the quantum Monte-Carlo (QMC) [l|, the exact di- 
agonalization Q , and the other methods such as the dy- 
namical cluster approximation (DC A) [1,01 have clarified 
many aspects of the cuprate superconductors, the prob- 
lems on the low energy excitations such as the non Fermi 
liquid behavior and the pseudogap state remain unre- 
solved theoretically because of the limited range of the 
inter-site correlations and the limitation of the momen- 
tum and energy resolutions in the numerical calculations. 

In order to describe long-range intersite correlations 
with high resolution in energy and momentum space, 
we have recently proposed a full self-consistent projec- 
tion operator method (FSCPM) for single-particle excita- 
tions [5| ■ The theory is based on the projection operator 
technique for the retarded Green function, and makes use 
of an off-diagonal effective medium to describe the non- 
local correlations efficiently. The off-diagonal matrix ele- 
ments of the effective medium are determined to be con- 
sistent with those of the self-energy of the Green function. 
We demonstrated on the basis of the half-filled Hubbard 
model on the simple cubic lattice that the long-range an- 
tiferromagnetic correlations cause shadow bands in the 
low-energy region and sub-peaks of the Mott-Hubbard 
bands in the strong interaction regime. 

In this proceedings, we apply the theory to the two- 
dimensional Hubbard model which is the simplest model 
of the cuprates, and investigate nonlocal effects on the 
momentum-dependent single-particle excitation spectra 
at zero temperature. We will demonstrate that the 



FSCPM describes well low-energy excitations of the 
doped two-dimensional Hubbard model with high reso- 
lutions and that the nonlocal correlations can cause the 
shadow bands as well as the shift of the Mott-Hubbard 
band splitting in the strong correlation regime. In partic- 
ular, we suggest that a Fermi surface instability should 
take place at the doping concentration 6% = 0.123. This 
might correspond to the 1/8 instability of cuprates. 

In the following section, we briefly review the FSCPM 
theory of single-particle excitations. In the theory we 
start from the Laplace transform of the retarded Green 
function. Introducing an off-diagonal effective medium, 
we expand the memoryf unction with use of the incremen- 
tal cluster expansion [6]. A remarkable point is that the 
off-diagonal effective medium is self-consistently taken 
into account up to the infinity so that high resolution 
is achieved in both energy and momentum. In Sec. Ill, 
we present numerical results of momentum dependent ex- 
citation spectra of the doped Hubbard model. In the last 
section we summarize a conclusion of the present work. 



II. SELF-CONSISTENT PROJECTION 
OPERATOR APPROACH 

We consider the retarded Green function, and adopt 
the Hubbard Hamiltonian H with nearest-neighbor 
transfer integral t and intra-atomic Coulomb interac- 
tion U. In the projection operator method 0, the re- 
tarded Green function Gij a (z) is expressed with use of 
the Laplace transformation as follows. 
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Here the inner product between two operators A and B 
is defined by (A\B) = {[A+,B]+), ( ) ([ , ]+) being 
the thermal average (the anti-commutator). a\ a {ai„) is 
the creation (annihilation) operator for an electron with 
spin a on site i. Furthermore z denotes the complex 



1 



energy variable z = uj + iS with 8 being an infinitesimal 
positive number, and L is a Liouville operator defined 
by LA — [H,A\- for an operator A. Here [ , ]_ is the 
commutator. 

The Green function is expressed as follows according 
to the Dyson equation. 



G l3 Az) = [(z-H»-A{z))-'] 
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and (A(z))ij a = Aij a (z) = U Gij a (z) is the self-energy 
matrix. Reduced memory function Gij a (z) is given by 



G ija (z) = Ul 
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The operator A\ a is defined by A\ a — a\ a 8ni^ a with 
8ni- a = n%-a — {ili— and (nj CT ) is the average electron 
number on site i for spin a. L is a Liouville operator 
acting on the space orthogonal to the space {laj^)}; L = 
QLQ, Q = 1 - P, and P = |a| CT ) (a| ff |. 

In the FSCPM theory, we introduce an energy- 
dependent Liouville operator L(z) for an effective Hamil- 
tonian with an off-diagonal medium Y,ij a (z); Hq(z) = 
Hq + ^ijcr( z ) a \<j a j<j- Here H is the Hartree-Fock 
Hamiltonian. The Green function Fij a {z) for the Liou- 
villcan L(z) is expressed as 



F ija (z) = [(z-H -t(z))- x } ija , (4) 

where {Y,(z))ij a = Y>ij a (z). It should be noted that the 
Green function F i j rj {z) becomes identical with Gij a {z) 
when 



(5) 



In order to obtain an explicit expression of the self- 
energy Aij a (z), we separate the Liouvillean L into L(z) 
and the remaining interaction part L\(z), i.e., L = L(z) + 
Li(z), and expand the memory function with respect 
to the interaction by using the incremental method @. 
In the present calculations, we take into account the in- 
tersite correlations within the pair-site approximation; 
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G iia (z) = GZ(z) + J2 AG uJ^) 

Gij a (z) 
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These terms 



Here AG^(z) 
are calculated from cluster memory functions defined 

by = (A\ a \{z-L {c \z))^A) a ) (c = 

L {c \z) = QL {c \z)Q and L^(z) is the Liouvillean 
for a cluster c embedded in the off-diagonal medium 

{t,lma(z)}. 

In the calculation of the cluster memory functions men- 
tioned above, we applied the form obtained by the renor- 
malized perturbation theory Q. 
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Here (Lj = U(l - 2{n i ^ a ))/xia, and Xia = 

(n i _ .)(l — (rij_ ff )). In the simplest approximation 

( G o Jija ( z ) is given by 
{G{^)ij a (z) = 
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Here (-Ho)ijcr is the Hartree-Fock Hamiltonian matrix jjere A- 



[Xia/( n i-o-}c(l - (ni_ CT ) c )]<5y + 1 — 5i- 



X (e,e',e") = (1 - /(e))(l - f(e'))f(e") + f(e)f(e')(l 
f{e")), /(e) is the Fermi distribution function, (nj CT ) c is 
the electron number for a cavity state defined by (n ia ) c = 



fdep%(e)f(e), and p™(e) = —ir~ l Im[(F c (z)~ 
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The coherent cluster Green function 

and 



U<T - F ija (z) is given by Eq. (gj), 

(S v ' (z))ijcr = J^ijcr(z) for sites (i,j) belonging to the 
cluster c. Note that a "cluster" c = (ij) does not mean 
that sites are nearest neighbors. Instead they may 
be far apart. 

Since we have truncated the higher-order terms in 
the expansions ([6]) and 0, the self-energy Aij a (z) — 
U 2 Gij a (z) depends on the medium T,ij a (z). We deter- 
mine the medium self-consistently from condition ([3]). 
Note that the present theory reduces to the projection 
operator CPA in infinite dimensions which is equiva- 
lent to the dynamical mean-field theory [Iol - [l3 |. 

The momentum dependent excitation spectra are cal- 
culated from the Green function 



G ka (z) 
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Aka(z) 

Here efc CT is the Hartree-Fock one-electron energy eigen 
value, and the momentum-dependent self-energy is calcu- 
lated via Fourier transform of the off-diagonal self-energy 
as Aka(z) = J2j A-jOa(z) exp(ik ■ Rj). Note that the the- 
ory yields the spectra with high resolution in both energy 
and momentum because it is based on the retarded Green 
function and the off-diagonal effective medium is taken 
into account up to infinity in distance. In the DCA, 
the spectra with high resolution are not obtained at low 
temperatures because it relies on the numerical analytic 
continuation and a small size of cluster embedded in an 
effective medium. 



III. NUMERICAL RESULTS 

We have performed the self-consistent calculations 
of excitation spectra for the two-dimensional Hubbard 
model on the square lattice in the nonmagnetic state 
at zero temperature. In the calculations of the memory 
functions (j^]), we adopted the Laplace transform, and cal- 
culated off-diagonal self-energies up to the 50th nearest 
neighbors self-consistently. 

Figure 1 shows calculated single-particle densities of 
states (DOS) for U = 8 in unit of |t| = 1 at half- 
filling. Note that the Mermin- Wagner theorem [TJJ tells 
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FIG. 1: The single-particle DOS at half-filling in the FSCPM 
(solid curve) , the SCPM-0 (dashed curve) , and the SSA (dot- 
ted curve). 



FIG. 2: Contour map of the momentum-dependent excita- 
tions along the high-symmetry lines at 8h = 0.14. Closed cir- 
cles with error bars are the QMC results at T = 0.33. Dashed 
curve shows the Hartree-Fock quasiparticle dispersion. 



us the antiferromagnetic state with the Neel tempera- 
ture Tn = K at half-filling, while we assumed here the 
nonmagnetic state. Thus the results presented here are 
not directly applicable at half-filling, but express the ba- 
sic behavior of excitations in the vicinity of half-filling. 
The DOS is compared with those of the single-site ap- 
proximation (SSA) and the diagonal self-consistent ap- 
proximation (SCPM-0) in which the off-diagonal effec- 
tive medium has been neglected, but the off-diagonal self- 
energies Aijff(z) have been taken into account. The result 
indicates that the van Hove singularity remains near the 
Fermi level in the vicinity of half-filling. We find that 
nonlocal correlations increase the splitting between the 
upper and lower Hubbard bands as compared with the 
SSA. On the other hand, the peaks of the both Hubbard 
bands are suppressed as compared with those obtained by 
the SSA and the SCPM-0. The increment of the splitting 
is explained by strong antiferromagnetic correlations. In 
fact the latters cause the splitting U + 2z^\J\ instead of 
U in the strong correlation regime. Here znn is the num- 
ber of the nearest neighbors, and J is the super-exchange 
interaction J = — 4|i| 2 /f/. The formula yields the split- 
ting 12, and explains the splitting in Fig. 1. 

Nonlocal correlations also cause small peaks at energy 
~ 3.0. This is interpreted as shadow bands of type 
e| DW (fc) = ±\/e| + ^ 2 d ue t° long-range antiferromag- 
netic correlations. Here £k is a quasiparticle band, A is 
an exchange splitting. 

When holes are doped, the spectral weight of the lower 
Hubbard band moves to the quasiparticle band as well as 
the upper Hubbard band. Furthermore the upper Hub- 
bard band shifts to the higher energy region. Figure 2 
shows an example of the momentum dependent excita- 
tion spectra in the optimum doped region. The excita- 
tions of the lower Hubbard band are damped, and those 
of the upper Hubbard band are enhanced. These re- 
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FIG. 3: Excitation spectrum at the Fermi energy for 
5h = 0.05. The dotted curve shows the Fermi surface in the 
Hartree-Fock approximation. 

suits of spectra are consistent with those of the QMC [l[ . 
Especially, the quasiparticle band near the Fermi level 
shows a good agreement with that of the QMC. 

Excitations at energy ut — determine the Fermi sur- 
face. We present in Figs. 3 and 4 the spectra at zero 
energy in the underdoped region (Sh — 0.05) as well 
as the overdoped region (8h = 0.20). Due to a strong 
damping of the lower Hubbard band, a hole-like Fermi 
surface appears in the underdoped region, while in the 
overdoped region an electron-like Fermi surface appears. 
The results are in good agreement with those based on 
the QMC combined with the DC A ||. 

Although the results of the FSCPM are consistent 
with those expected from the QMC+DCA at finite tern- 
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FIG. 4: Excitation spectrum at the Fermi energy for 
&h = 0.20. The dotted curve shows the Fermi surface in the 
Hartree-Fock approximation. 
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FIG. 5: The density of states at 8h = 0.123 where the quasi- 
particle peak is just on the Fermi level. 



peratures, there are a few points showing the discrep- 
ancy. The marginal Fermi liquid behavior is one of 
them. The marginal Fermi liquid is characterized by the 
imaginary part of the self-energy being proportional to 
max(|wL T) Q. In the recent 4x4 QMC+DCA calcula- 
tions [la ] , the marginal Fermi liquid behaviors are found 
at finite temperatures. The result for U — 6 suggests that 
the marginal Fermi liquid remains even at zero temper- 
ature at doping concentration St — 0.15, where the peak 
of the single-particle DOS is just on the Fermi level. This 
is rather close to the result of the SSA, 6^ = 0.153 in our 
calculation {i.e., the projection operator CPA). 

We point out that the peak at St is connected to that 
in Fig. 1 at zero doping. With the hole doping, the peak 
with van Hove singularity sinks first below the Fermi level 



because of the transfer of the spectral weight from the 
lower Hubbard band to the upper Hubbard band. But, 
it gradually approach to the Fermi level with further dop- 
ing, and the peak is again located on the Fermi level at 

In the present calculations of the FSCPM, the critical 
doping concentration at which the quasiparticle peak is 
located just on the Fermi level is 6% = 0.123 for U = 8 
as shown in Fig. 5. Furthermore electrons are the Fermi 
liquid instead of the marginal Fermi liquid at St because 
we find that the imaginary part of the momentum depen- 
dent self-energy in the quasiparticle energy region is pro- 
portional to lu 2 instead of It is remarkable that the 
critical concentration — 0.123 is close to the 1/8 stripe 
instability found in the experiments [l6l Il7j . Though the 
present results indicate the Fermi liquid behavior, it is 
possible that the Fermi surface causes an instability at 
5^ because of the peak just on the Fermi level, so that the 
normal state may change to the stripe phase at S^ = 1/8. 



IV. CONCLUSION 

We have investigated the nonlocal effects of elec- 
tron correlations on the excitation spectra of the two- 
dimensional Hubbard model at zero temperature on 
the basis of the full self-consistent projection operator 
method (FSCPM). The FSCPM takes into account the 
long-range intersite correlations in a self-consistent way 
making use of the off-diagonal effective medium, and 
yields the spectra with high energy and momentum res- 
olution. In the strong Coulomb interaction regime (i.e., 
the case of U — 8), we found that the long-range in- 
tersite antiferromagnetic correlations create the shadow 
band excitations. The strong antiferromagnetic correla- 
tions also enhance the Mott-Hubbard splitting by about 
2znn|</| where znn (J) is the number of the nearest 
neighbors (the super-exchange interactions energy). By 
comparing with the QMC and the QMC+DCA results, 
we verified that the FSCPM quantitatively describe the 
momentum dependent excitation spectra in the doped re- 
gion. We also found that the nonlocal correlations shift 
the critical concentration 5^ at which the quasiparticle 
peak is located just on the Fermi level from 6^ = 0.153 
(the SSA) to <5* = 0.123, while the QMC+DCA at fi- 
nite temperatures suggest S^ — 0.15 at zero temperature. 
Furthermore the electrons at 6% — 0.123 obtained by the 
FSCPM are the Fermi liquid, though the QMC+DCA 
suggested the marginal Fermi liquid even at zero temper- 
ature. It is plausible that the Fermi surface instability 
takes place because of the singular peak on the Fermi 
level. We speculate that such an anomaly corresponds to 
the 1/8 stripe instability and that it may cause the non 
Fermi liquid behaviors at finite temperatures as found in 
the QMC+DCA in the vicinity of S* h . Details of the DOS 
and Sfr as a function of the Coulomb interaction U will 
be published elsewhere. Theoretical calculations for pos- 
sible instabilities at are left for future investigations. 
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